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Abstract 

The binding energy of H to a (10,0) carbon nanotube is calculated 
at 24, 50, and 100% coverage. Several different bonding configura- 
tions are considered for the 50% coverage case. Using the ONIOM 
approach, the average C-H bond energy for the most stable 50% 
coverage and for the 100% coverage are 57.3 and 38.6 kcal/mol, 
respectively. Considering the size of the bond energy of H 2 , these 
values suggest that it will be difficult to achieve 100% atomic H 
coverage on a (10,0) nanotube. 

I. INTRODUCTION 

We recently 1 computed the binding energy of hydrogen atoms to the side walls 
of a (10,0) carbon nanotube using the ONIOM method 2-4 . The first C-H bond was 
21.6 kcal/mol. The average C-H bond strength for the first two hydrogen atoms 
was 40.6 kcal/mol and for the first four hydrogens 47.9 kcal/mol. While there is 
an increase in the bond strength with number of hydrogens, the values were still 
small compared with the H-H bond at the same level of theory (109.8 kcal/mol). 
Therefore we suggested that even at high hydrogen coverage the C-H bonding would 
be endothermic or only slightly exothermic. Considering the interest 5 7 in storing H 
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using carbon nanotubes, in this manuscript we report on the results of calculations 
that model high H coverages on the same (10,0) carbon nanotube. 

II. MODEL AND METHODS 

The initial coordinates of the 20 A segment of the (10,0) nanotube were generated 
using the code of Han 8 . The dangling bonds at the ends were tied off with hydrogen 
atoms, yielding a C 200 H 20 tube. For the 100% coverage, one additional H atom is 
bound to each carbon atom, thus yielding a C 200 H 220 species, which is shown in 
Fig. 1. For 24% and 50% coverages three random coverage patterns were considered. 
In each of these, the H atoms were added in pairs. The first H atom was added to 
one of the bare carbon atoms in a random manner. The second H atom was then 
added to one of the bare first nearest neighbors in a random manner. If there were no 
first nearest neighbors, the first H placement was rejected. This approach was used 
to break individual C-C if bonds, and hence maximize the remaining if bonding in 
the tube. 

The observation that F on the outside 9 and I on the inside 10 of carbon nanotubes 
have distinct patterns, rather than random coverages, leads us to investigate severed 
higher symmetry structures for the 50% H coverage; these cases are shown in Figs. 2 
to 6. Fig. 2 shows the case where the H atoms are parallel to the axis of the tube 
and evenly spaced around the tube. In Fig. 3, the hydrogens are also along the axis 
of the tube, but the pairs of lines of H atoms are adjacent. These two structures 
are denoted “lines” and “pairs of lines”, respectively. In Fig. 4, the hydrogens atoms 
spiral around the tube; this structure is denoted as “spiral”. In Figs. 5 and 6, there 
are rings of H atoms; in the first, denoted “rings”, the hydrogens are evenly spaced, 
while in the second case, the rings of hydrogen atoms are in pairs, which is denoted 
as “pairs of rings”. 

The AMI and two level ONIOM approach 2-4 are used. The ONIOM approach 
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combines the Universal Force Field 11 (UFF), for the low-level treatment, with the 
hybrid 12 B3LYP 13 functional, for the high-level description. The 4-31G basis set 14 
is used in conjunction with the B3LYP calculations. The AMI calculations are per- 
formed using a modified version of Gaussian 94, where damping is used to obtain 
convergence. ONIOM calculations are performed using Gaussian 98 15 . 

The geometries are fully optimized at the AMI and ONIOM levels of theory. In 
the ONIOM calculations, twenty-four carbon atoms, at the center of the nanotube, are 
used for the high-level treatment. The boundary hydrogen atoms and the chemisorbed 
hydrogen atoms are also in the high level treatment. The atoms included in the high 
level treatment are the same as used in our previous treatment 1 of the low coverage 
case. We should note that more recent work has shown 16 that capping the carbon 
nanotubes instead of terminating the dangling bonds with hydrogen atoms appears 
to dramatically speed up the geometry optimization process. 

As noted in our previous study, we encounted problems with local mirfima for the 
UFF description of the carbon nanotube, and therefore we use the B3LYP energies 
instead of the ONIOM energies, since the B3LYP energies were insensitive to the 
UFF solution. That is, we only use the molecular mechanics approach to constrain 
the shape of the high-level fragment. 

III. RESULTS AND DISCUSSION 

We first consider the results obtained at the AMI level of theory — see Table I. 
The average H binding energy for three 24% random coverages span a range of 
about 5 kcal/mol. Since there are 48 C-H bonds, the stability of the tubes vary 
by 243 kcal/mol. The average H binding energy for the three random 50% coverages 
spans a somewhat smaller range, and the values are similar to those obtained for 
the 24% coverage. We can conclude that the AMI binding energies for the 24 and 
50% coverages are similar, but it would require many more runs to more accurately 
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determine the average bond energies for these two cases and how it changes between 
them. However, as discussed below, we can get some insight into the bonding using 
only these few cases. 

In light of patterns observed for F and I with nanotubes, rather than consider 
more random coverages, we studied five higher symmetry patterns, which are shown in 
Figs. 2-6 and their average H binding energies are summarized in Table I. Three of the 
higher symmetry structures are more stable than the random coverages. Considering 
that there are 100 C-H bonds, the change in the average H binding energy means 
that the three higher symmetry structures are several 100 kcal/mol more stable than 
the random structures. 

The lower stability of the random structures arises for two reasons: 1) there are 
areas without H atoms that contain an odd number of carbons, which dramatically 
degrades the 7 r bonding, and 2) to form a good C-H bond, the carbon must sp 3 
hybridize, which results in the carbon bulging out of the tube. For the random tubes, 
there are areas of very high H coverage, and clearly the tube cannot deform such that 
every atom can form a good sp 3 hybrid bond with H. 

For the higher symmetry cases, all regions devoid of H atoms contain an even 
number of carbon atoms, so the 7r bonding is not degraded for this reason. For the 
pairs of lines tube, two rows of carbons bulge out of the tube allowing the formation 
of good C-H bonds. This deformation of the carbon tube leaves the bare rows of 
carbon atoms nearly planar; thus this configuration has both good C-H bonds and 
good C-C 7r bonds. In the spiral, the deformation leaves a spiral of good C-H bonds 
and a spiral of C-C re bonds. In the lines configuration, the geometry appears very 
favorable for good bonding, but the C-C 7r bonds are isolated, which removes the tt 
conjugation, and hence this is less favorable than the pairs of lines or spiral. In the 
two ring forms, the curvature of tube weakens the C-C tt bonding, thus these forms 
are unfavorable, even though they distort to form good C-H bonds. 
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We studied the two most stable 50% coverage structures, namely the spiral and 
pairs of lines, using the ONIOM approach. We computed the average H binding 
energy of the 12 hydrogens included in the high-level treatment — see Table I. While 
we do not compute the H binding energy of the hydrogens in the low-level treatment, 
their effect on the geometry of the tube is included. The pair of lines ONIOM value 
is very similar to the AMI result, while for the spiral, the ONIOM value is somewhat 
smaller than the AMI result. On the basis of our low coverage work, part of the 
difference between the AMI and ONIOM values can probably be attributed to the 
fact that many of the H atoms in the high level treatment are near the boundary 
between the high and low-levels of theory. However, it seems unlikely that this can 
account for the entire difference, and some, if not most, of the difference must be 
attributed to limitations in the methods used, with the AMI treatment having a 
larger uncertainty than the ONIOM approach. 

Using the ONIOM approach, we computed the average H binding energy for 100% 
coverage, and we find that this value is significantly smaller than the 50% coverages. 
For the 100% coverage, the carbon atoms still form a good tube structure since any 
deformation that improves one C-H bond will weaken another. This means that for 
the 100% coverage the carbons cannot change their hybridization to enhance the C-H 
bond, which results in a much weaker C-H bond than the 50% coverage, where half 
of the C atoms can bulge out of the tube to maximize the C-H bonding. 

Despite any limitations in the methods used, it is clear that the formation of 100% 
coverage will be a very endothermic process (remember that the H-H bond energy is 
109.8 kcal/mol at this level of theory). The computed binding energies for the 50% 
coverage cases are sufficiently close to one half the H-H bond energy that it might 
be possible to achieve this level of coverage in an exothermic process starting with 
H 2 and nanotubes. The formation of significantly higher than 50% coverages in an 
exothermic process seems unlikely since higher coverages would require a deformation 
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of the carbons in the tube that would weaking some of the existing C-H bond and 
hence result in a smellier average H binding energy. 

It should be noted that eill of our discussion is for the case of H atoms on the 
outside of the tube. It might be possible to achieve high coverage by bonding some 
H atoms on the inside of the tube and some on the outside. For example, one half 
the C atoms bulge out and bond to H atoms and the other half of the C atoms bulge 
in and bond to H atoms on the inside of the tube. However, this approach faces the 
problem of opening the tubes and getting H atom inside. 

IV. CONCLUSIONS 

We use the AMI level of theory to obtain some insight into the bonding of H atoms 
to a (10,0) carbon nanotube. The carbon atoms that bond to the H atoms bulge out 
of the tube to improve the sp 3 hybridization, which enhances the C-H bond. For 50% 
coverage we have found two structures that have half of the C atoms bulged out of the 
tube to form good C-H bonds, while the remaining carbons atoms can still have good 
7 r bonding. Random coverages are found to be much less favorable because not all 
atoms have a favorable orientation for either C-H or C-C x bonding. Random coverage 
also tends to lead to bare regions with an odd number of carbon atoms which is not 
ideal for x bonding. Using the ONIOM approach, the average C-H bond energy was 
computed for the two most favorable 50% coverage cases as well as for 100% coverage. 
It is possible that the formation of 50% coverage is thermodynamically favorable, but 
it is concluded that the formation of 100% coverage will be very endothermic. 
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TABLES 


TABLE I. Summary of hydrogen binding energies, in kcal/mol. 


System 

Average H binding energy 

AMI 

24% random case 1 

53.80 

24% random case 2 

52.08 

24% random case 3 

48.74 

50% random case 1 

51.11 

50% random case 2 

50.67 

50% random case 3 

48.91 

50% pairs of lines 

58.29 

50% spiral 

57.62 

50% lines 

56.11 

50% rings 

51.51 

50% pairs of rings 

48.54 

ONIOM 

1 H a 

21.6 

2 H’s° 

40.6 

4 H’s a 

47.9 

50% pairs of lines 

57.29 

50% spiral 

52.64 

100% 

38.65 


a Taken from reference 1. 
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FIGURES 

FIG. 1. The C 200 H 220 tube used to model 100% hydrogen coverage on a (10,0) carbon 
nanotube. 

FIG. 2. The C 200 H 120 tube used to model the “lines” form of the 50% hydrogen coverage 
on a (10,0) carbon nanotube. 

FIG. 3. The C 200 H 120 tube used to model the “pairs of lines” form of the 50% hydrogen 
coverage on a (10,0) carbon nanotube. 

FIG. 4. The C 200 H 120 tube used to model the “spiral” form of the 50% hydrogen cov- 
erage on a (10,0) carbon nanotube. 

FIG. 5. The C 200 H 120 tube used to model the “rings” form of the 50% hydrogen coverage 
on a (10,0) carbon nanotube. 

FIG. 6. The C 200 H 120 tube used to model the “pairs of rings” form of the 50% hydrogen 
coverage on a (10,0) carbon nanotube. 
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